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INTRODUCTION 


This paper shows the results of an in-house and contractual 
effort to better define the parameters associated with thermal contact 
conductance. Data for contact conductance vs. applied pressure, and 
the corresponding graphs are shown for samples of 304 Stainless Steel, 
AZ31 Magnesium, 6061-T6 Aluminum and Copper. 

For a more thorough discussion of the contractual work, the 
reader is referred to the final report of this contract (NAS8-5207). A 
paper covering the results of this contract will also be presented at 
the AIAA 1st Annual Meeting and Technical Display June 29 -July 1964, 
at the Sheraton Park Hotel, Washington, D. C. 

In addition to the contractual interface data, an attempt is made 
to define the observed change of slope of 6061 -T6 and 2024-T4 Aluminum 
when the data are plotted on log -log graph paper. It is shown that by 
deforming cones, hemispheres, and ellipses, a similar change of slope 
occurs. It is concluded that these models might possibly represent 
"scale-up" replicas of the macroscopic points of contact of two mating 
aluminum surfaces. 

A reference list is included which is a revision and extension 
of the bibliography the author handed out at the February meeting. It 
contains many previously unknown Russian references. 



EXPERIMENTAL PROGRAM 


A study of the problems in the early stages of the thermal con- 
tact conductance work, has indicated a need for experiments designed 
to (1) aid in the understanding of the heat transfer mechanism, (2) pro- 
vide data to verify existing ana'yses, (3) provide data to aid in the de- 
velopment of new analytical methods. 

Subsequently, a thermal contact conductance apparatus suitable 
for use in vacuum was developed which would permit accurate measure- 
ment of thermal conductance as a function of contact pressure. As 
opposed to the flat plate apparatus used in the investigations reported 
by Fried, the principal investigator of this study, this apparatus utilized 
cylindrical columns to minimize flatness deviations under load. 

Thermal Test Apparatus 

A schematic of the test apparatus is shown in Fig. 1. Figure 2 
shows the heat flow section of the apparatus, with a specimen in place, 
without the radiation shield. 

The samples consisted of two metallic cylinders having a diameter 
of 5. 08 cm (2 in. ), and a length of 7. 62 cm (3 in. ) each. Each sample 
was instrumented with four copper constantan thermocouples to determine 
the axial temperature gradient due to the uniform heat flux passing be- 
tween the electric heater and the liquid-cooled sink. 

Contact pressure could be varied by means of a stainless steel 
bellows, pressurized in accordance with the desired load. The load was 
measured using a strain gage load washer on the heat sink side (Fig. 1). 
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The entire assembly was installed in a bell jar vacuum system 
with a right angle cold trap, utilizing a 4 -inch oil diffusion pump pre- 
ceded by a roughing pump to achieve a vacuum of 10 4 mm Hg (1. 33 
x 10 * newton/m 2 ) or better. 

The heat source utilized in this test was a 100 -watt electric 
resistance element embedded in the main heater assembly which is 
gaurded by a ring heater and a rear guard heater, as shown in Fig. 1. 
This system is arranged such that there exists no temperature difference 
between the main heater and the guards. Each is separately controlled, 
so that all thermal energy from the main heater has only one direction 
to go — into the test sample. In order to monitor this system, thermo- 
couples were fastened to the several surfaces seeing each other. 

Minimum cross-sectional area supports, made of tubes (Fig. 1), 
were used between the rear guard and the main heater, in order to 
minimize heat leak errors, even though the facing surfaces were kept 
at the same temperature. The desired range of temperature differences 
between potential heat lead points were kept at AT's of 1 °C or less 
in order not to exceed l/2 of 1% heat flow errors. Initially, these 
temperature differences were controlled by use of a deviation amplifier , 
but experience indicated that manual control, with proper judgment, 
resulted in less time delay between steady -state points. 

The allowable temperature differences were dictated by the 
amount of heat passing through the test sample, since high heat fluxes 
through the sample permitted higher heat losses from the heater while 
permitting the percentage losses to remain the same. 
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The heat flux was determined by measuring the regulated d-c 
power input (i. e. , voltage and current), using precision instruments. 

In addition to this, the hot heater resistance was obtained by momen- 
tarily turning off the power. In order to eliminate leadline losses in 
the calculation, the ratio of heater winding resistance to total system 
resistance was measured and a correction applied to all readings. An 
ESI bridge having an accuracy of + 0. 05% was used. 

A check was performed on the adequacy of the heat flow measure- 
ment by determining the thermal conductivity of a piece of ARMCO 
iron. The measured value came within 2% of the nominal value which, 
considering all possible variables, is quite good. If we were to per- 
form only thermal conductivity measurements, this accuracy could 
probably be improved. However, for conductance measurements, with 
their many sources of error, the cost of improving this system is not 
quite worth the effort at present. 


Temperature Measurement 

Considerable attention was paid to accurate temperature 
measurement techniques in order to minimize possible measurement 
errors, since the quality of the temperature measurement directly 
affected the quality of the interface thermal contact conductance obtained. 
Thermocouple junctions were made of 40 -gauge copper -con stantan pre- 
cision grade thermocouple wire. This grade of wire has a nominal 
tolerance of + 0, 3 °C over the range of interest, but has been found by 
experience to be considerably better. Junctions were made by mercury 
pool arc welding techniques. 
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The thermocouples were installed in the test samples in 2. 54- 
cm deep holes, to place the junction at the cylinder axis. The junction 
was embedded with Eccobond 56C, an epoxy base cement having a 
thermal conductivity equal to that of stainless steel. In order to assure 
that the thermocouple bead actually contacted the sample at the cylinder 
centerline, a 0. 33 -cm diameter hole was drilled at the desired axial 
thermocouple location and a tube of the same material as the sample 
was inserted with the thermocouple installed. This method had the 
advantage that there was less likelihood of drill runout when the hole 
was drilled. It also permitted more positive installation and location 
of the thermocouple junction. The only exception to the matching of 
material was that an aluminum tube was used with the magnesium 
sample. This was not expected to result in an error because: (1) the 

thermocouple junction was in contact with the sample magnesium, and 
(2) the thermal effect of different material was not adverse because of 
the higher thermal conductivity of the aluminum. This would not result 
in a delay to reach thermal equilibrium. 

The choice of 40 -gauge thermocouple wire was dictated by the 
desire to minimize conduction losses. Experience with several hundred 
thermocouples from such wire (purchased from Thermo-Electric Co.) 
with no adverse emf characteristics led to the selection of this diameter. 
The question as to the proper response of the thermocouples when 
embedded in the samples in a vacuum was circumvented by use of the 
Eccobond 56C, a fairly free -flowing epoxy cement inserted and packed 
around the thermocouple bead and wire. Thus, the bead was hermetically 
isolated from the surrounding atmosphere. 



To assure proper response of these thermocouples, they were 
placed in a constant temperature oven after being installed in the 
sample and the consistency of the temperature readings was checked. 
Out of over 60 thermocouples tested, only 4 were found to require 
corrections in the computation of conductances for the range of temper- 
atures of interest (25-50°C). Particular attention was paid to the 
precision with which the axial distances between thermocouples were 
controlled, since the axial distance vs. temperature plots were used 
to project the temperature gradients to the interface and thus obtain 
the interface temperature difference. 

The constriction resistance effects at and near the interface 
require that thermocouples be located in the undisturbed region in 
order to correctly project the temperature gradient. Since only the 
sample half interfaces are of interest, the heat source and heat sink 
interfaces with the samples had high vacuum silicone grease applied as 
a heat transfer promoting device. Thus, no significant constriction 
effects resulted at these interfaces. 

The temperature difference, AT , is based on the temperature 
obtained experimentally, which are then extrapolated to the interface. 
The accuracy with which this AT can be obtained is a function of the 
accuracy with which the temperature gradient in the sample can be 
obtained. For high values of contact conductances the AT usually 
was quite low. Conversely, for low value s of conductance the AT 
was high. Since a high AT resulted in a higher percent accuracy, 
the relative percent accuracy of contact conductance obtained was 
constant. A representative temperature gradient curve is shown in 



Fig. 3. Of the thermocouples used in the samples, each had its own 
cold junction. Their emf was read on a Leeds and Northrop K-3 
potentiometer, with individual couples switched by means of a transfer 
switch. Figure 4 shows the vacuum system, thermocouple recorder, 
power supply, and instrument panel. 

Surface Finish Measurements 

One significant area of interest, which strongly affects the 
thermal contact resistance is the surface finish of the interface. 

Surface finish, by definition, can include surface roughness as well 
as waviness, which is described by Clausing and "microscopic and 
macroscopic effects, " and by F enech as "primary and secondary 
wa vine s s . " 

In addition to the small asperities which constitute the rough- 
ness, a machines surface can have larger peaks and valleys which 
constitute the waviness. The direction parallel to the ridges and 
valleys of the waviness is called the lay direction. 

A Taylor -Hobson "Talysurf" stylus type profilometer was used 
to obtain single-line profiles of the various surface finishes prepared 
to this program. Due to difficulties of operating ari in-house "Talysurf" 
instrument, all but one pair of samples (Nos. 15 and 16) were inspected 
after thermal contact conductance tests were completed. Thus, any 
deformation of asperities, which may have taken place during tests 
would, therefore, be observable. However, it is not very likely that 
any such effects could be observed, because the "Talysurf" trace is 
merely the record of a stylus motion following the contours of the 
surface in a straight line. 
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Any asperity, deformed or otherwise, on either side of this 
straight line would, therefore, not be recorded. Although there is no 
certainty that a trace parallel to or in continuation to an existing trace 
will resemble the existing trace, there will be a similarity of char- 
acteristics, provided the character of the surface is taken into con- 
sideration. For example, in the case of machined surfaces, traces 
should be taken in the direction of tool motion as well as in the perpen- 
dicular direction. Particular attention should be paid to lathe -turned 
finishes at the profile through the center of the surface, because of 
the non-flatness of the surface at that point. Figure 5 shows typical 
"Talysurf" traces through the center of a machined surface for a 
copper sample. 

The traces as shown, do not represent a true pictorial repre- 
sentation of the surface, because of the scale differences. These 
asperities appear to be much more severe than they are in reality. 
Nevertheless, the traces do provide a significant amount of useful 
information and provide an excellent means for comparison of surface 
finishes. 

As a result of the length of the stylus travel (1 . 27 -cm max. ) 
which is adjustable, and the use of the optical flat attachment, flatness 
deviations can also be observed. This is due to the fact that the stylus 
motion, relative to an optical flat, is recorded. 

An additional feature of the "Talysurf" profilometer is its 
ability to provide a centerline average (CLA) roughness reading, by 
means of an electronic integrator circuit, for any surface of certain 
minimum length. Centerline average (CLA) is also known as arithmetic 
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average (A A) and runs somewhat lower than the corresponding root- 
mean-square (RMS) reading. The latter gives more weight to the 
larger deviations from the centerline. 

Flatness measurements were made using a surface plate and 
a dial indicator reading 2. 5 micrometers, (0. 001 inches) which per- 
mitted estimation of half divisions (1. 3 micrometers). The dial indi- 
cator point was set at the sample center and the dial was set at zero. 

With the dial indicator fixed, the sample was moved so that the point 
traveled to the interface edge, reading the vertical deviation at the 
center, one -fourth diameter and at the edge. 

This was done at mutually perpendicular diameters. A secondary 
check was made initially by holding the sample fixed and moving the 
dial indicator support stand. No significant differences were observed 
between the two methods. Plus readings indicated high spots, whereas 
minus readings indicated low spots. Results are shown in Table I in 
which the maximum values are presented. It should be noted that these 
values are the maximum from a fictitious plane, i. e. , the datum plane 
as described in the next major section, "Deformation Experiments. " 
Thus, there may occur some matching of interfaces having deviations, ' 
which could result in a test assembly of better mating than would be 
expected on the basis of individual reading. For example, samples 
3 and 4 could have a cumulative flatness deviation of only +1,2x10 6 
meters if they fitted into each other. 
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Thermal Test Results 


The material and the important surface properties of the test 
samples are shown in Table I. These include roughness, Rockwell 
hardness, flatness deviation and type of surface preparation. Actual 
data for these surfaces are shown in Table II. 

Stainless Steel 304 

Figure 6 shows the results of the stainless steel interface tests. 

Of interest is the large difference in conductance at the maximum con- 
tact pressure. The flatness deviation of the 0. 30 micrometer (RMS) 
roughness samples was 1. 3 micrometer, whereas the 1.2 micrometer 
(RMS) roughness sample had a flatness deviation of approximately 1 . 5 
micrometer, at best, and 3.8, at worst, depending on surface matching. 

Of interest is the curvature of the fine finish contact conductance 
curve whose behavior was confirmed by the descending load curve. 

Hysteresis could be observed for this specimen for the loading -unloading 
cycle. 

In contract, the coarse finish sample curve shows no hysteresis 
and is almost linear . 

It is of particular interest to note and compare these two curves 
in Fig, 6 with the corresponding results of Clausing. The resemblance 
of the Clausing results with Stainless Steel 303, for approximately the 
same degree of flatness deviation, to our results is remarkable. The 
importance of the approximate similarity of flatness deviation, as 
opposed to a marked difference in roughness (Clausing, 3 micro -in for 
both versus our 12 and 50 micro -in) is demonstrated well in this experiment. 


10 



Magnesium 


Figure 7 shows the results for Magnesium AZ31B, a widely- 
used magnesium alloy. These samples, which had lathe -turned inter- 
faces exhibited a rather unusual reversal of expected performance. 

The coarse finished surfaces exhibited higher thermal contact conduc- 
tances than did the fine finished interfaces. One possible explanation 
would be the greater effect of a surface film on a fine -finished surface 
versus that on a coarse -finished surface. Oxide films and tarnish were 
visible on both sets of samples, since two months had elapsed between 
machining and use. The reason for conjecture that a film will have a 
lesser effect on a coarse surface finish, is that the fewer sharper 
ridges of this finish will result in higher loads per unit area and cause 
the film to break. Another, and perhaps more plausible, reason is 
the relatively large flatness deviation for both sample pair, but that 
the sample assembly may have resulted in a greater mismatch for the 
poorer performance. 

It is of interest to note that Clausing obtained higher conductances 
for similar material having lower values of flatness deviation and much 
lower surface roughness. 

Aluminum 

The resultant conductance versus pressure curves are shown in 
Fig. 8. It is of interest to note that there was ho significant difference 
in the values of contact conductance for the two surface finishes con- 
sidered. The results for the finer (0. 3 -micrometer RMS) finish 
6061-T6 Aluminum should have been higher than for the coarse (1.4- 
micrometer RMS) finish, since the former had lower values of flatness 
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deviation. At present, no explanation can be found for this behavior. 

The general shape of this curve conforms to that shown by Clausing 
for 2024 Aluminum, with the conductance somewhat lower at maximum 
pressure . 

Copper 

A test for electrical grade copper (OFHC oxygen-free, high- 
conductivity copper) was performed, because the only available data 
(Jacobs and Starr) indicated linear variation of conductance with load 
at moderate loads, whereas, most other materials change in a non- 
linear manner in that pressure region. As can be seen in Fig. 9, the 
curve is not linear at low pressures, but does appear to be linear at 
higher contact pressures. It is also of interest to note that no hysteresis 
could be observed for this copper joint. 

General Remarks 

The results for specific metal joints are discussed under their 
respective headings. This section discusses common-ground observa- 
tions. 

When conductance versus pressure is plotted on log -log paper, a 
Curve (as shown in Figs. 11-13) results, which is somewhat different 
from earlier observed and expected results. Initially, a slope of one - 
half to two -thirds was expected for elastic behavior as discussed in 
another section of this paper. However, plots of data obtained in this 
study indicate a definite two -regime behavior with a pronounced point 
of change in slope. The exact reason for this change in slope has not 
yet been defined, except to show that it possibly represents the change 
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from purely elastic to elastic -plastic deformation behavior. This is 
discussed in the next section dealing with an experimental study of 
this phenomenon. 


DEFORMATION EXPERIMENTS 

The three models (2024-T4 Aluminum) described in this paper 
are shown in Fig. 10. The cone and hemisphere models were 2. 54 cm 
(1 in. ) in diameter and 1. 27 cm (0. 5 in. ) in height. The ellipse semi- 
major axis was 1 . 27 cm with its semi -minor axis being . 950 cm (.375 in. ). 

The models in Fig. 10 (column 1) were placed between two 
flat plates of a steel press with a piece of pressure -sensitive paper 
placed on their tops and bottoms and a load Pj was applied. A typical 
piece of the pressure -sensitive paper appears below the models. The 
blackened area is the deformed area for that particular load. After 
each specified loading, another paper was placed on the model. Over 
the entire range of loading from 0-250, 000 Newtons (0-60, 000 pounds), 
the deformed area remained circular, as indicated by the blackened 
area on the paper, and the deformed model. The diameter of this 
blackened area was measured several times and an average taken, 
thus leading to the recorded deformed area data in Table III. The tests 

O 

were performed at room temperature (293 K), 

The height of the model was measured by a dial micrometer 
placed between the two steel plates. The models in columns 2, 3,4, 
and 5 of Fig. 10 were subjected to specific loads, and the areas and 
heights were compared to those of the previously described models in 
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which the load was cycled. No appreciable difference was noticed and, 
thus, the cycling of loads had produced little work hardening of the 
models. 


As soon as the data were plotted, it was observed that an 
interesting resemblance existed between the published thermal interface 
data and the deformation of the model. Of particular interest is that of 
the area/height deformation versus loading when compared to the thermal 
interface conductance as a function of its mechanical loading. Figure 
11 shows data of the models compared on a log -log plot with that of 
Fried and of Clausing. In an attempt to bring the data into the same 
order of magnitude, the following expression was used; 



P=P. 

x 



(1) 


where 


k 



conductivity of the models 

deformed area of the model at load (P^) 

height of the model at load (P^) 

computed conductance of the models 


to compute a representative thermal conductance. It must be strongly 
emphasized that the plotted data in Figs. 11-13 taken from the Fried 
and Clausing reports should not be used in computation. This data 
has been shifted in magnitude for better visual observation. 

It is particularly interesting that both the interface data and 
the model data experience a change of slope at certain loading values. 
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The factor that appears to cause this change of slope in the model data 
is the dependence of the deformed area on the loading. This became 
evident when the area versus the loading was plotted. The contribution 
of the model height versus loading did not undergo this sudden change. 

This critical point of loading at which the slope changes shall, hereafter, 
be designated P c j for the interface data and P c m for the model data. 

As can be seen from Fig. 11, the values of P c j and P c j^ do 
not coincide. This might be partly explained by a temperature depen- 
dence. In comparing P c m with Pci °f Clausing, it is to be noted 
that the models were at 293 °K (70° F) while Clausing reported mean 
interface temperatures of approximately 386°K (234° F) for eight 
interfaces. Figure 12 is a plot of the data reported in this paper for 
6061 aluminum and the computed model data. This mean interface 
temperature (T-^) was approximately 301° K (82 °F), this value being 
the average of all the Tj and T 2 values of interface. Since for this 
sample T m was near that of the model temperature, it appears that 
the slope change at P c j is nearer the value of that of the models P c ]y[ 
than the corresponding Clausing data. However, this comparison is 
not totally valid since the metals are different. This leads to the question 
of whether P c j is dependent on the mean interface temperature. If 
P c j is attributed to the changes of the physical properties of the metal, 
it would appear reasonable that its value should be lower for higher 
mean interface temperatures. Thus, it would appear that 

K d .at (P c l) <** ( T~ ) ( 2 ) 

x m 
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If all the load values of the deformation models are divided by 
the corresponding deformed area, pressure values are recorded which 
are consistently near the yield strength of the metal, as can be expected 
for permanent deformation. 

It seems that there are other factors which influence P c j 
for the Clausing data. If the eight data groups are plotted, then P c j 
appears at different load values for each specimen. This is partly 
shown by two curves of Fig. 13 . 

When all the eight samples of 2024-T4 Aluminum values are 
averaged and plotted, Fig. 13 shows that the two -slope regime is 
again evident. As can be seen, this corresponds to the included data 
for the models. 

In order to study the functional relationship of the curves a 
computer program for best fitting the data to an equation was formed. 
This equation corresponds to the form presented earlier and is 

Q 

h = A + BP . The data from Clausing, data reported in this paper, 
and the deformation model data show similar values of the exponent 
c both before and after the change of slope. The values of A/Y, h, 

A, B, do not coincide because the data used for the best fit curve were 
of different units as reported in the respective reports. On Figs. 12 
and 13, only the functional notation has been shown for comparison. 

The best fit curves are: 

1 . Model data 

^ = - 32. 90 + 0. 57 P°' 72 
from P = 0 to 10,000 pounds. 
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A -a ] S4 

— = - 37.47 + 3.03 x 10 4 P 

from P = 10, 000 to 60, 000 pounds, 
where 

A = deformed area in inches 2 
Y = deformed height in inches 
P = load in pounds 


2 , 


6 06 1 -T 6 Aluminum (Fig. 12) 

h = 9. 73 x 10 +2 + l. 17 x 10 +3 

from P= 1 0 . 2 to 41 9 p. s , i . 

h * 1. 14 x 10 +2 + 7. 00 x 10“ 2 

f r om P = 419 to 1, 117.0 p.s.i, 
where h is given in BTU/hr ft 2 


p0.09 


P 


1 .61 


O 


F . 


3. Average data of Clausing for eight samples 
of 2024-T4 Aluminum 

0 91 

h = 35.41 + 7. 59 p 
from P = 1 0. 4 to 67. 0 p. s. i. 

h = 1 68 . 1 + 2. 14 p 1 ' 16 
from P = 67. 0 to 986. 0 p. s. i. 
where h is given in BTU/hr ft 2 °F. 


CONCLUSIONS 

1. The importance of the flatness deviation effects on thermal joint 
conductance has been demonstrated. 
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2. The proposed models, based on the elastic deformation relations 
of Hertz appear to provide an approach to understanding the heat 
transfer mechanism. This is represented by the approaches of 
Clausing and this paper. 

3. Better surface definition methods are required. 

4. More experimental data of suitable accuracy is needed to arrive at 
(a) semi -relations and (b) statistical correlation. 
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TABLE III 


DEFORMED AREA AND HEIGHT OF 
MODELS AS A FUNCTION OF APPLIED LOAD 


Model 

Load 

Area 


Height 

kilonewtons 

kilopounds 

meters 2 
x 10' fo 

Inches 

millimeters inches 

All Models 

0 

0 

0 

0 

12.700 

.500 

Cone 

, . 445 

. 100 

. 107 

. 002 

12.421 

. 489 

Cone 

1.335 

.300 

.324 

. 005 

12. 294 

. 484 

Cone 

2.224 

.500 

• 636 

.010 

12.065 

.475 

Ellipse 

2.224 

.1500 

1. 140 

.018 

12. 598 

.496 

Hemisphere 

2. 224 

.500 

1. 265 

. 020 

12. 624 

.497 

Cone 

3. 559 

.800 

1. 140 

. 018 

11.938 

.470 

Cone 

5. 338 

1.200 

1.534 

. 024 

11.735 

.462 

Cone 

6.672 

1. 500 

2.027 

.031 

11. ‘582 

. 456 

Ellipse 

6.672 

1.500 

2.634 

. 041 

12. 497 

.492 

Cone 

8.896 

2.000 

2. 588 

. 040 

11,481 

452 

Hemisphere 

8.896 

2.00 

3.426 

. 053 

12.497 

.492 

Ellipse 

11. 121 

2. 500 

3.973 

.062 

12. 370 

. 487 

Cone 

13.345 

3.00 

3.694 

.057 

11.024 

.434 

Ellipse 

15.569 

3. 500 

5.451 

.084 

12.319 

.485 

Cone 

17.793 

4. 000 

4. 560 

.071 

10. 693 

. 421 

Hemisphere 

17.793 

4.000 

6.936 

. 108 

12. 319 

.485 

Ellipse 

22.241 

5.000 

7.946 

. 123 

12. 090 

.476 

Hemisphere 

26.689 

6.000 

9.813 

. 152 

12. 167 

.479 

Cone 

31. 138 

7.000 

7. 240 

. 112 

10. 033 

. 395 

Ellipse 

33.362 

7. 500 

12. 067 

. 187 

11.862 

.467 
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Model 

Load 

Area 

Height 


kilonewtons 

kilopounds 

meter sz 
x 10 -5 ' 

inches 

* millimeters 

inches 

Cone 

44. 482 

10.000 

9.810 

.152 

9.601 

.378 

Ellipse 

44, 482 

10.000 

14.234 

.221 

11.557 

.455 

Hemisphere 

44. 482 

10.000 

15. 329 

.238 

11.887 

.468 

Ellipse 

55. 603 

12. 500 

18. 241 

.283 

11.252 

.443 

Cone 

66.723 

15. 000 

15. 328 

.238 

8.560 

. 337 

Ellipse 

66.723 

15.000 

23,430 

. 363 

10.922 

.430 

Cone 

88.964 

20. 000 

13 . 155 

.369 

7. 595 

.299 

Elipse 

88.964 

20. 000 

28. 199 

.437 

10. 135 

. 399 

Hemishere 

88.964 

20.000 

31. 142 

... 483 

11. 125 

.438 

Cone 

111.206 

25.000 

34.071 

.528 

6. 731 

.265 

Ellipse 

111.206 

25.000 

37.826 

. 586 

9. 347 

.368 

Hemisphere 

111.206 

25.000 

40.677 

. 631 

10. 643 

.419 

Cone 

133.447 

30.000 

42. 888 

.665 

6. 020 

.237 

Ellipse 

133.447 

30.000 

47. 480 

.736 

8.458 

.333 

Hemisphere 

133.447 

30.1000 

48.071 

.745 

10. 109 

.398 

Cone 

155.688 

35.000 

54. 806 

.849 

5. 385 

.212 

Ellipse 

155.688 

35. 000 

59.981 

.930 

7.595 

.299 

Hemisphere 

155.688 

35. 000 

58.013 

.899 

9. 550 

. 376 

Ellipse 

177.928 

40. 000 

70.077 

1.086 

6. 756 

.266 

Hemisphere 

177.928 

40. 000 

67.477 

1. 046 

8.992 

.354 

Hemisphere 

200. 170 

45. 000 

78.413 

1. 215 

8.433 

. 332 

Hemisphere 

2221411 

50. 000 

90.174 

1.400 

7. 976 

.314 

Hemisphere 

266. 893 

60.000 

123. 948 

1.921 

T. 163 

.282 




FIGURE I. VARIABLE PRESSURE -THERMAL INTERFACE 
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Fig. 5 


THERMAL CONTACT CONDUCTANCE VS. APPLIED PRESSURE 
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THERMAL CONTACT CONDUCTANE VS. APPLIED PRESSURE 
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THERMAL CONTACT CONDUCTANCE VS. APPLIED PRESSURE 




Fig. 10 
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Fig. 13 P ( MECHANICAI LOAD IN KILONEWTONS) 
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